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Abstract 

We discuss the propagation of neutrino wave packets in a Lense-Thirring space- 
time using a gravitational phase approach. We show that the neutrino oscillation 
length is altered by gravitational corrections and that neutrinos are subject to he- 
licity flip induced by stellar rotation. For the case of a rapidly rotating neutron star, 
we show that absolute neutrino masses can be derived, in principle, from rotational 
contributions to the mass- induced energy shift, without recourse to mass generation 
models presently discussed in the literature. 
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Recent experimental evidence from Superkamiokande [1] and SNO [2] has 
significantly endorsed the claim [3] that neutrinos undergo flavour oscillations 
in vacuum due to the difference of their rest masses. This discovery, however, 
makes no pronouncements about the intrinsic physical properties of neutrinos, 
particularly whether they exist as plane waves or wave packets. In addition, 
the actual calculation of the oscillation length for solar neutrinos assumes a flat 
space-time background. While this assumption seems reasonable for neutrino 
propagation in our solar system, it may be unwarranted when the gravitational 
source is a neutron star or a black hole. 
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The purpose of this work is to investigate the contribution of gravitation to the 
flavour oscillations of a two-neutrino system in a Lense-Thirring space-time 
background [4], assuming a wave packet description of the mass eigenstates 
in momentum space. This approach [5], so far untested in gravitational prob- 
lems, offers an intuitively satisfying description of neutrinos. It relates directly 
the physical characteristics of the gravitational source, such as its stellar tem- 
perature and density, to the size of the wave packet. It also lends itself well 
to sensing distortions in space-time because of the common and intrinsically 
non-local nature of both curvature and wave packets. The choice of Lense- 
Thirring space-time is particularly relevant for our purpose, since rotational 
effects may induce a helicity transition while the neutrino is in transit. 

A particularly interesting feature of our approach is the introduction of the 
gravitational phase that leads to a direct spin-gravity coupling interaction 
within the neutrino's wavefunction. It was shown [6,7,8,9] that the gravita- 
tional phase 
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enters the description of quantum particles in external gravitational fields in 
a way that is essential and independent of their spin. In (1), 7^ = — i]^ is 
the metric deviation, and p M and L a @ are momentum and orbital angular mo- 
mentum operators of the free particle. Application of (1) to closed space-time 
integration paths gives rise to a covariant Berry's phase [10] with consequences 
for particle interferometry [11]. Equation (1) also yields the particle deflection 
predicted by general relativity in the geometrical optics approximation [12]. 
When applied to fermions, <3>g and the spin connection T M reproduce all those 
gravitational-inertial effects that have been either observed [13,14] or pre- 
dicted [15,16], and predicts, in particular, the non-conservation of helicity [17] 
for strictly massless fermions. 

To show how $g acts on a neutrino propagating in vacuum, we first start 
with the covariant Dirac equation [i^(x) (V M + iT^) — m/h] ip(x) = 0, where 
m is the neutrino rest mass, and V M is the usual covariant derivative. We 
use geometric units of G — c — 1 [18], so that all physical quantities can be 
described in units of length. Then the corresponding Dirac Hamiltonian in 
Lense-Thirring space-time is 
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where r = yfx 2 + y 2 + z 2 , M and i? are the mass and radius of the gravi- 
tational source, Q is its angular velocity, and L z = xp y — yp x is the orbital 
angular momentum operator in the z-direction. The gravitational phase can 
be introduced by means of the transformation if)(x) — * exp (i^q/H) ip(x) and 
the new Hamiltonian takes the form H = H + H$ G , where 



H^ G = <x ■ (V$ G ) + (V*$ G ) (3) 

is a first-order correction. We then treat (3) as a perturbation for a two mass- 
species neutrino system. While the approach taken is based on standard quan- 
tum mechanics, in order to apply (1) it is necessary to assume that the average 
of all possible paths reduces to the integration path of $ G . Because some terms 
of H are non-diagonal with respect to spin, the wave packet approach should 
also make known whether the neutrino mass eigenstates are subject to a he- 
licity transition in the gravitational field of a rotating source. 

Adopting the Dirac representation, the matrix element is (ip(r)\H^ G \ip(r)), 
and we assume for the wave packet description 



Wr)> = (2^5/2 / rf3 ^( fc ) e ife - p |tf(*)> ± , ( 4 ) 

where £(fc) = €(k x )€(k y )£(k z ) is the normalized Gaussian wavefunction of 
width <7 P and centroid (fe), and 
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The normalized four-spinor \U(k))± is 



(5) 
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where E = y/p 2 + m 2 is the energy associated with the spinor. By symmetriz- 
ing over the exchange between k and k', the explicit construction of the matrix 
element shows contributions due to both a spin diagonal term and a helicity 
transition, in the form 
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where (V^cls = (V 2 ) [(V m $g) (fe) + (V^g) is the symmetrized form 
of the gravitational phase gradients, and V = (Att/3) (r 3 — R 3 ) is the volume 
of spatial integration from the star's surface. It is clear from the last line of 
(7) that there would be no helicity transition contribution if we considered 
strictly plane waves, since a non-zero transition amplitude requires different 
momentum components within the wave packet to interact with the Pauli spin 
matrices to yield a non-zero result. Choosing the spin quantization axis to be 
along the neutrinos' direction of propagation, the helicity transition element 



is 



(=F|cr|±) = [cos# cosy? =F i simp] x + [cos 6* simp ± i cosip] y — sin#z, (8) 

where the upper sign refers to the transition from negative to positive helicity. 

The explicit calculation of the integrals of (7) in spherical coordinates is pro- 
hibitively complicated by the coupling of the oscillatory functions to the ampli- 
tudes. However, an approximate, but analytic expression for (7) can be found 
by integrating over both position and momentum space after performing a 
2nd-order Taylor series expansion with respect to to and an llth-order Taylor 
expansion of the oscillatory functions. It can then be shown that (7) has the 
form 
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where to = m/(p) = m/{h{k)) and q = (k)/a p , with Cj(q,r) and Dj(q,r) 
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Fig. 1. Functions which determine Eq. (9) for 10 MeV neutrinos, given a r = 10 
kpc neutron star source of M = 1.5 AT© and R = 10 km. 

as dimensionless functions that can be tabulated. The first line of (9) is the 
contribution due to the diagonal components of (ip(r)\H<p a \ijj(r)) , while the 
second line is due to the off-diagonal components and refers to the helicity flip 
contribution of the perturbation. 

For this paper, the gravitational source under consideration is a rapidly rotat- 
ing 1.5M neutron star with R = 10 km and = 1 kHz, at a distance of r = 10 
kpc. Figure 1 contains a list of plots for the functions described in the matrix 
element (9). For this choice of parameters, it is clear from a comparison of Fig- 
ures l(a)-l(c) and Figures l(d)-l(f) that \C 2 (q,r)\ < \d(q,r)\ < \C (q,r)\ 
and \D 2 (q,r)\ <C \Di(q,r)\ <C \D (q, r)| for all choices of q. This suggests the 
trend towards convergence of the series due to expansion with respect to m. To 
justify the truncation of the Taylor expansion of (7), we note from analyzing 
a one-dimensional analogue of the problem that the Gaussian functions which 
are present in the integrand have the effect of damping out the contribution 
of the higher-order expansion terms beyond the 11th order. A plot of this se- 
ries expansion matches precisely with that due to the corresponding product 
function of the exact sinusoidal function with a one-dimensional Gaussian, 
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Fig. 2. Corresponding functions for Eq. (12), assuming the same source and neutrino 
energy conditions. 



except in the tail regions of the Gaussian envelope. In those isolated regions 
of parameter space, we estimate an error of 10-15% in the amplitude and are 
confident that this degree of error has no significant bearing on our results. 

We stress some interesting features of (9). The most obvious one is the presence 
of terms linear in m, due to the Taylor expansion of (ijj(r)\H<s> G \ijj(r)) for both 
the spin diagonal and off-diagonal terms. This fact has interesting implications, 
as shown below, in the calculation of the energy shifts. Given (8), we know that 
the helicity transition contribution to (9) is due entirely to the z-component 
of the transition amplitude, since the x- and ^/-components average out to 
zero over all spatial angles. Furthermore, the non-zero contribution is due 
entirely to the rotation of the source, which induces the helicity transition. 
Because of the sin 9 term, propagation of a neutrino along the iz-axis suggests 
that its initially prepared helicity state remains fixed throughout its motion. 
This is a sensible result that is consistent with our present understanding of 
spin-rotation coupling, since the rotational term in the Lense-Thirring metric 
tends to act like a uniform magnetic field that orients the particle spin either 
parallel or antiparallel to the field strength. This result also agrees with the 
fact that, barring quantum fluctuations about the classical integration path of 
(1), V^q = in the Lense-Thirring field, which implies helicity conservation 
[17]. 

To determine the mass-induced energy shift for the neutrino oscillation length, 
we note that the off-diagonal elements of (9) contribute to a second-order 
effect in time-independent perturbation theory, where the unperturbed energy 
eigenvalues Eq come from H ip{x) = Eq ip(x) for Eq « J (p) 2 + m 2 — — (p) + 
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By virtue of (10), we can calculate the energy shift for a given neutrino, leading 
to a new energy of E± m Eq + (AE')t , where 



(A^=(±|^ G |±)± l( Jf^ + _ ) 
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This leads [19,20] to the final expression for the neutrino oscillation length 
L osc . = 2tt/ (e^-E^), where 
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and 
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The plots of (14) and (15) are listed in Figure 2. 
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It is clear from Figures 1 and 2 that the functions become extremely large 
for q — > 0, corresponding to very large momentum spread in the wave packet 
for a given neutrino energy, which then rapidly decay to zero for large q. 
In particular, the plots show that the helicity transition terms in (12) will 
dominate as q — > 0, which is consistent with (7), since a large momentum 
spread in the wave packet is required to sense the differential rotational effects 
within a localized region of space-time near the star's surface. 

For solar neutrinos and those due to supernovae [20] , the expected wave packet 
widths in momentum space are ha p ~ 10 -5 MeV and 2 x 10~ 2 MeV, respec- 
tively. The corresponding values for q 10 6 and q ~ 50 indicate that grav- 
itational corrections have negligible effect on the neutrino oscillation lengths 
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for these scenarios. However, the situation is quite different when applied to 
rapidly rotating neutron stars. For the neutron star parameters we consider, 
M/r « 7.175 x lCT 18 and MVtR 2 /r 2 « 7.751 x 1CT 36 . In order to predict 
a 1% correction to the value of Am^ = m 2 , — m\ as determined by solar 
neutrino experiments, we require \F 2 (q,r)\ 1.3667 x 10 33 , which suggests a 
choice of q ~ 2.1 x 10 -5 from Figure 2(b), and implies a wave packet width 
of ha p ~ 4.762 x 10 5 MeV. This prediction corresponds very well to the cal- 
culated width of ha p ~ 3.260 x 10 5 MeV for neutrinos emitted from neutron 
stars, as determined by a mean-free-path calculation [20], assuming a stellar 
temperature of 3 x 10 6 K [21], and an effective stellar density of 10 11 g/cm 3 
averaged over the neutron star's expected core and surface densities. Our re- 
sults therefore show that helicity flip likely plays a role in the case of neutrinos 
propagating in the field of rotating neutron stars. 

From Figure 2(a), it also follows that \Fi(q, r)| 1.2150 x 10 38 for the same 
choice of q as when applied to Figure 2(b), and so the contribution in (12) 
due to linear mass difference Am 21 = m 2 — m\ is not negligible. This result 
suggests the possibility of performing, in principle, a parameter fit of (12) 
for suitable choices of q, Am 21 , and Am 2 i so that we can infer the absolute 
neutrino masses entirely from observations, and without any reference to the 
mass generation mechanisms presently under consideration in the literature 
[19]. In practice, however, such an undertaking would require an enormously 
large neutrino flux and large counting rates to obtain statistically significant 
measurements. To our knowledge, there are no experimental facilities available 
that can meet these requirements. Nonetheless, the theoretical possibility of 
determining absolute neutrino masses by this technique makes for an interest- 
ing consideration in the development of future neutrino observatories. 
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